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Abstract: A nonperturbative, purely numerical, solution of the reduced 
Rayleigh equation for the scattering of p- and s-polarized light from a 
dielectric film with a two-dimensional randomly rough surface deposited 
on a planar metallic substrate, has been carried out. It is found that satellite 
peaks are present in the angular dependence of the elements of the mean 
differential reflection coefficient in addition to an enhanced backscattering 
peak. This result resolves a conflict between the results of earher approxi- 
mate flieoretical studies of scattering from this system. 
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In the earliest analytic |[T] and computer simulation ||2][3) studies of the multiple scattering of 
light from clean one-dimensional randomly rough surfaces of perfect conductors or of penetra- 
ble media, the focus was on the phenomenon of enhanced backscattering. This is the presence 
of a well-defined peak in the retroreflection direction in the angular dependence of the intensity 
of the light that has been scattered incoherently (diffusely). 

In subsequent work on the multiple scattering of light from free-standing or supported films 
with a one-dimensional randomly rough surface that support two or more guided waves, new 
effects were discovered |4 |. These include enhanced transmission, which is the presence of a 
well-defined peak in the anti-specular direction in the angular dependence of the intensity of 
the light transmitted through the film |5|. Perhaps more interesting was the discovery of satel- 
lite peaks in the angular dependence of the intensity of the light scattered from or transmitted 
through the film. These are well-defined peaks present on both sides of the enhanced backscat- 
tering and enhanced transmission peaks, respectively, that arise from the coherent interference 
of guided waves with the frequency of the incident light, but with different wavenumbers ||6). 

It should be noted, however, that the prediction of these satellite peaks was first made in the 
context of the scattering of electromagnetic waves from a dielectric film containing a random 
distribution of volume scatterers |7|, rather than from a randomly rough surface, when the 
thickness of the film is small compared to the mean free path of the electromagnetic wave in 
the random medium. 

In analytic [8 -10] and computer simulation calculations fTTJ|T3) of the multiple scattering 
of light from clean two-dimensional randomly rough surfaces of perfect conductors and pen- 
etrable media, enhanced backscattering was observed in the results. However, when attention 
turned to the scattering of light from a perfectly conducting surface coated with a dielectric 
film, conflicting results were obtained. In these studies the dielectric -perfect conductor inter- 
face was assumed to be planar, while the vacuum-dielectric interface was assumed to be a 
two-dimensional randomly rough interface. In the first of these studies Kawanishi et al. | [T4) 



applied the stochastic functional approach to this problem and found no evidence for satellite 
peaks in their results. They suggested that the ensemble averaging of the intensity of the scat- 
tered field restores isotropy in the mean scattering plane, and thereby eliminates the occurrence 
of special scattering angles at which satellite peaks could occur In subsequent work in which 
the reduced Rayleigh equation for scattering from this structure | T0][T5) was solved in the form 



of expansions of the amplitudes of the p- and s-polarized components of the scattered field in 
powers of the surface profile function through terms of third order, satellite peaks were found. 
However, the contribution to the scattering amplitudes associated with the third-order term was 
larger than that from the first-order term for the roughness and experimental parameters as- 
sumed in that work. It is therefore possible that these values fell outside the ranges for which a 
perturbative solution of the reduced Rayleigh equation is reliable. 

Although satellite peaks were observed in experiments carried out by Mendez et al. 1 16] that 
utilized the double passage of polarized light through a random phase screen, the experimental 
conditions were sufficiently different from those studied theoretically in Refs. 1 10 T4|, that 



these results could not be used to support the predictions of either of these studies. 

In an effort to resolve the issue of whether satellite peaks do or do not exist in the scattering 
of light from a rough dielectric film deposited on the planar surface of a metal, in this paper we 
carry out a nonperturbative, purely numerical, solution of the reduced Rayleigh equation p7) 
for the scattering of p- and s-polarized light from a structure consisting of a dielectric film 
deposited on a metal substrate when the dielectric-metal interface is planar, while the vacuum- 
dielectric interface is a two-dimensional randomly rough interface. This is an approach that 
was used successfully in recent calculations of the scattering of p- and s-polarized light from 
a two-dimensional randomly rough interface between a dielectric and a metal |[T7][T8|, which 
prompts its application to the present problem. 

The system we study consists of vacuum (ei) in the region X3 > d + C, (x||), where X|| ~ 
{x\,X2,Q)', a dielectric film (£2) in the region < X3 < d + tl, (x||); and a lossy metal (£3) in 
the region X3 < 0. The surface profile function 1^ (x||) is assumed to be a single-valued function 
of X|| that is differentiable with respect to xi and X2, and constitutes a zero-mean, stationary, 
isotropic, Gaussian random process defined by 

c(x||)c(x;i)U52w 



(1) 

The angle brackets here denote an average over the ensemble of realizations of the surface 
profile function, and 5 = (1^^ (x||))'''^ is the rms height of the surface roughness. 

The electric field in the vacuum [x3 > d+ll, (x||)] is the sum of an incident field and a scat- 
tered field, E(x;f) = [E(x| 'B)mc +E(x|fi))sc] exp (— iof), where 



E(x|a))i„c = I ^ [k|| ai {k\\ ) +t?,k\\\ Bp (k|| ) + (xj x k|| ) B, (k,, ) 
X exp (ik|| • X|| — iai )x3) 



(2a) 



E(x|to)sc =y (^1^ [q||«i(?||)-X3^||]^p(q||) + (x3 xq||)As(q||)j ^^^^ 

X exp (iq|| • X|| + iai (^11 )x3) , 

while the subscripts p and s denote the p-polarized and s-polarized components of these fields 
with respect to the local planes of incidence and scattering. A caret over a vector indicates that 
it is a unit vector, and the vector k|| is defined as k|| = {ki ,k2,0) (with similar definitions for qii 



and p||). The functions Cli{qn ) (/ = 1 , 2, 3) are defined by 



CO 



1/2 



Re a,- (^11 ) > 0, Ima, (^||) > 0. 



(3) 



A linear relation exists between the amplitudes Aa(q||) and B^(k||) {a,j3 = p,s), which we 
write as 



Aa(q||)=L^a/3(q|||k||)B/3(k|| 



(4) 



where R^p is the scattering amplitude for incident /3 -polarized light scattered into a-polarized 
light. The convention we use with respect to the polarization subscripts is 



R(qii|k||) 



%(q|||k||) /?ps{q|||k,|) 
7?sp(q|||k||) /?ss(q|||k||) 



(5) 



It has been shown by Soubret et al. pO) and Leskova p9) that the scattering amplitudes 
\j^ap (q|| |k|| )] satisfy the matrix integral equation 



/ 



(27r)- 



-M(p,||q[|)R(q|||k||) = -N(p|||k||), 



(6) 



called a reduced Rayleigh equation because it is an equation for only the scattered field in the 
medium of incidence, and not for the fields in the film and in the substrate. The effects of the 
latter two fields are contained in the elements of the matrices M(p|| |q|| ) and N(p|| |k|| ). With the 
shorthand notation a(<7|| , co) = ci{qH ), the elements of these matrices in the forms obtained by 



Leskova | |19j are 

A^pp(piilqii) = [pii^ii + «2(/?||)(P|| ■q\\)aiin'>] 

^(«2(P||)-«i(^||)|p||-q||) 



: rp(/:>|| ) exp (-i [ai {p\\)- a\ iq\\ )] d 
-[/'||^||-«2(p||)(P||-q||)o:i(?||)] 
: Ap(p||)exp(i [a2ip\\) + ai (^||)] d 



[Pll^ll -«2(p||)(P||-q||)o:i(?||)] 

[«2(/5||) + o:i(?||)] |P||-q|| 



^ps(P|||q||) = --«2(P||)(P|| xqii 



a2ipii) + ai{q\\) 

3 



'r,(„).xp(-.K(„,,-a,(„)].)'>"''*-"'<'''l"-'") 



Ap(/?||)exp(i [a2(/J||) + ai(^||)] d 



^sp(P|||q||) = ^ (P|| xq||)3o:i(?|| 



^s{p\\ ) exp (-i [a2 (pii)- ai (^|| )] d) 



a2{p\\) + ai{q\\) 
/(a2(/?||)-ai(^l|)|p||-q||) 



As (pii ) exp (i [a2{p\\ ) + ai (^|| )] t/) 



a2{p\\) + ai (^11 ) 

A^ss(P|||q||)= "^(P||-q||) 



«2 



a2(p||) + «i(^ll 



(7a) 



a2{pii) - ai(qii) (7b) 
[a2(/'||) + ai(?||)] |P|| -q||)\ 



«2(/5||)-ai(^l|) (7c) 
[«2(p||) + ai(?||)] |P|| -q||)\ 



^ / -r / ^ / M A'f(«2(p||)--ai(?||)|P||-q||) 
r, (;,,, ) exp (-1 [a2(P|| ) - 0=1 (^11 )] d) ^;(^^^),^;(^;') (7d) 

,A/ ^ ^ ^, ^ ,1 A'^(-[«2(/?||) + o:i(?||)]|P||-q||)^ 

+ As(/:'ll)exp(i[a2(p||) + ai(?||)Jt/) — ' 



and 



^Pp(P|||k||) ^ - [/?||^|| -a2(/?||)(p|| •k||)ai(^||)] 

/(a2(/?||) + ai(^ll)|P|l -k||) 



: rp(/7|| ) exp (-i {ai (/?|| ) + ai {k\\ )] d) 
- [pii^ii +a2(p||)(p|| •k||)ai(^||)] 
: Ap (pii ) exp (i [a2 (/?|| ) - ai (^11 )] <i) 



a2(p||) + ai (^11 ) 

[pil^ll +a2(p||)(p|| •k||)ai (^||)] 

/(- [a2(/.||)-ai(^||)] |p||-k|| 



a2(p||)-ai(^ll) 



A^ps(p|||k||) = --a2(/?||) (P|| xk|| 



3 



V ^ ^ ( -l i ^ , (1 M A'^(«2(/?||) + ai(^:||)|P||-k 



Ap(;,||)exp(i[a2(;,||)-a,(fc||)]^) 



Afsp(P|||k||) = (p|| xk||)3ai(^||) 



«2 



A'ss(Piilkii)-^ P||-k 



and 



as well as 



(8a) 



a2(P||) + ai(^||) (8b) 
/(- [a2(p||)-ai(^||)]|p||-k||)\ 



(/,,, ) exp (-1 [a2 (/.„) + ai (^|| )] d) ^^^^^^ ^ '^^(^^^ (8c) 

Ki ^ (-1 i ^ n M ^ [«2(/'|| ) - «i (^||)] IP|| - k||) ^ 
As(;'||)exp(i[a2(;'||)-ai(^l|)Jt/)— ' 



ir f ^ ( ■\ I ^, M A'^(«2(/'||) + 0:i(^||)|P||-k||) 

X (^r.(p||)expH[a2(;,||) + a,(^||)]^) aM) + o^^{f^d ^'''^ 

,A/ ^ /T ^ ^ M^^^(-[«2(P||)-«i(^||)]|P||-k||)\ 
+ A. ) exp (i [a2 ) - ai (A:,, )] ^) ) ■ 

In writing Eqs. (|7]i and (|8]l we have introduced the functions 

rp(/?||) = e2a3(p||,t(j) + e3a2(p||,t(j) (9a) 
Ap(p||) = £2a3(p||,«)-e3a2(pi|,«) (9b) 



T,{p\\ ) = ajipw , a) + a2(p|| , a) (10a) 
As(/?||) = a3{p\\,(o)-a2{ppco), (10b) 

/ (rlQll) - / d\ exp (-iQ|| -xii) exp [-iyC (x,,)] . (11) 

The scattering amplitudes [Rap (q|| |k|| )] play a central role in the present theory because the 
mean differential reflection coefficient, an experimentally measurable function, can be ex- 
pressed in terms of these amplitudes. The differential reflection coefficient (dR/dCls) is defined 



such that {dR/dQ.s)dQ.s is the fraction of the total time-averaged flux incident on the surface 
that is scattered into the element of solid angle dD.^ about the scattering direction {9s, ^i). Since 
we are studying the scattering of light from a randomly rough surface, it is the average of this 
function over the ensemble of realizations of the surface profile function that we need to calcu- 
late. The contribution to the mean differential reflection coefficient from the incoherent (diffuse) 
component of the scattered light, when incident light of )3 polarization whose wave vector has 
the projection k|| on the mean scattering surface is scattered into light of a polarization whose 
wave vector has the projection q|| on the mean scattering surface, denoted ((9/?„^/(9ni).^^^j^, is 
given by 



(12a) 
(12b) 
(12c) 
(12d) 



do.,, 1 


^ incoh 


1 Ver»«?(?ii) 

~5 47r2 c ai(/t||) 






^ incoh 


1 e]!^ CO a^iq\\) 
^5 47r2 c ai{k\\) 




dRsp' 




1 1 CO a'^iq\\) - 


do.,, 1 


^ incoh 


S47z^y/e{ c ai{k\\) . 


dR,,' 


^ incoh 


1 ^iCoa^{q\\) 
S47z^ c ai(/t||) 





/?pp(q|||k||)|')-|(/?pp(q|||k||))| 
^ps(q|||k||)|')-|(/?ps(q|||k||))|' 

-K^sp(q|||k||)) 
(^ss(q|||k||)) 



^sp(q|||k||)| 



^ss(q|||k||) 



where S is the area of the plane X3 = covered by the rough surface. The two-dimensional 
wave vectors k|| and q|| are defined in terms of the polar and azimuthal angles of incidence 
(00, 0o) and scattering (0^,0^), respectively, by k|| — ^/e\{co/c) sin0o (cos0o,sin0o,O) and 
q|| — -yer(fi)/c) sin 0i (cos sin 0i,O). Thus these wave vectors in Eq. ([T2]i are restricted to the 
domains < s/e{{co/c) and <7|| < ^ei"((o/c) of the qxqi plane. 

Up to now Eq. (|6| has been solved by small-amplitude perturbation theory through terms 
of third order in the surface profile function ||T0 20 1. Here we present results for the mean 
differential reflection coefficient and for the full angular distribution of the intensity of the 
scattered light obtained by a nonperturbative, purely numerical solution of Eqs. (|6l)-([TT]l, as 
described in Ref. p7) . This was done by generating a realization of the surface profile func- 
tion numerically on a grid of A^^ points within a square region of the x\X2 plane of edge L, so 
that the (linear) sampling interval was Ax = L/Nx- A two-dimensional version of the filtering 



method used in 1 17 21 1 was used to generate the profile function ||T3|. The function /(7|Q||) 
was then evaluatMby expanding the integrand in powers of the surface profile function C(x||), 
and calculating the Fourier transform of CH'^n) ^y the fast Fourier transform. In evaluating the 

integral over q[| in Eq. (|6]l the infinite limits were replaced by finite ones: (ci\ +^2) ''^^ — G/2- 
The Nyquist sampling theorem requires that \qi \ and \q2\ be smaller than Qc = n/Ax |22 p. 
605]. The components of the vector p|| — q|| entering /(7|p|| — q||) lie in the interval [- 
so we have chosen Q = Qc- A grid with a grid constant A^i — Aq2 — Aq 
structed within the circular region of the qiqi plane where (g'j 



-Q.Ql 
In/L was con- 
< 2/2. The integral over 

this region in Eq. (|6]l was carried out by a two-dimensional version of the extended midpoint 
method |22., p. 161] and the values of R^^^ (q|| |k|| ) were calculated for values of q|| at the points 
of this grid for a given value of k||, which was also a point on this grid. The resulting matrix 
equations were solved by LU factorization and backsubstitution. The values of Rg.^ (q|| |k|| ) and 
(q|| |k|| ) P were then calculated for Np realizations of the surface profile function. An arith- 
metic average of the A^^ results for each of these functions yielded the averages (R-afi (q|| |k|| )) 
and (|^a|g(q|| |k|| )|^), from which the incoherent contribution to the mean differential reflection 
coefficients were calculated according to Eq. ( [T2] i. 



We apply this approach to the scattering of p- and s-polarized plane waves, whose wave- 
length is A = 633 nm, incident from vacuum (ei = 1) on a dielectric film (£2 = 2.6896 + 0.01i) 
coating a silver surface (£3 = — 18.28 + 0.481i) |23J. The mean thickness of the film is 
d = 0.756A = 478.5 nm. The roughness of the vacuum-dielectric interface is characterized 
by a two-dimensional version of the West-O'Donnell power spectrum p4) given by ||8| 

^(|k||l) = ^T^0(|k|||-^-)0(^+-|k|||), (13) 

where Q{x) is the Heaviside unit step function, and = 0.82(a)/c), k+ = 1.97((b/c). The 
rms height of the surface roughness was assumed to be 5 = X/AQ = 15.82 nm, the surface 
was discretized on a grid of resolution zVjci = Ax2 = 0.123A = 77.6 nm and the edge of the 
(quadratic) surface was L = 55A = 34.8 [a.m. 

The contribution to the mean differential reflection coefficient ((97?„p(q|| |k||)/(9ns). ^ 

from single-scattering processes [second order in C(x||)] is proportional to g{\i\\\ ~k|||) 18 1. 
Since the power spectrum ( pjj ) is identically zero for |k|| | < there is no contribution to the 
mean differential reflection coefficient from the light scattered incoherently by single-scattering 
processes when the wave vectors q|| and k|| satisfy the inequality |q|| — k|| | < The contri- 
bution to ((9/?c(jg(q|||k||)/(9i2j).^^^j^ when this condition is satisfied is due only to multiple- 
scattering processes, including the enhanced backscattering peak and the satellite peaks. These 
features are more clearly visible in this case because they do not ride on a large background 
due to single-scattering processes. This is the reason that the calculations whose results are 
presented here were carried out on the basis of the power spectrum ( pj) . 

In Fig. 1(a) we present the contribution to the mean differential reflection coefficient from 
the light scattered incoherently as functions of the polar scattering angle 9^ for the in-plane 
(0i ~ ^0 — 45°) co-(p— s-p, s^s) and cross-(p— >s, s->p) polarized scattering when a p- or s- 
polarized plane wave is incident on the dielectric surface at angles of incidence (0o,0o) given 
by (0.74°, 45°). (In figures showing in-plane or out-of-plane scattering, we depart from the com- 
monly accepted principle of not using negative polar angles, in that we allow for negative 0(.). 
An arithmetic average of results obtained for A^^ ~ 11, 165 realizations of the surface profile 
function was carried out to produce these figures. In Fig. 1(b) we present the analogous results 
for out-of-plane (0, — 0o±9O°) scattering when the roughness and experimental parameters 
have the values assumed in generating Fig. 1(a). 

In the results depicted in Fig. 1(a) [1(b)] single-scattering processes give no contributions to 
the mean differential reflection coefficient for -53.8° < 0, < 56.4° (-55.08° < 0, < 55.08°). 
In both figures a well-defined enhanced backscattering peak is seen in the retroreflection direc- 
tion. In addition, in Fig. 1(a) additional peaks are seen on both sides of the enhanced backscat- 
tering peak in the s s co-polarized scattering contribution to the mean differential reflection 
coefficient. These peaks are identified as satellite peaks. 

We base this identification on the following consideration. It was shown in Q that in the in- 
plane co-polarized scattering of light of frequency o from a one-dimensional randomly rough 
surface of a film system when the plane of incidence is perpendicular to the generators of the 
surface, satellite peaks occur at scattering angles given by 



c 



sin 0; ' = - sin 00 ± 1= [qm{(0) - q„{(0)] . (14) 



The wave numbers qi{co),q2{co), ■ ■ ■ -.qNico) are the wavenumbers of the guided waves sup- 
ported by the film structure at the frequency of the incident light. Not all of the peaks predicted 
by Eq. ( [T4| may be present in the mean differential reflection coefficient. This happens when the 



absolute value of the right-hand side of Eq. ( 14 1 is greater than unity. Then the corresponding 




-90 -60 -30 30 60 90 



Fig. 1. (a) The contributions to thie mean differential reflection coefficient as functions of 
the polar scattering angle ds from the in-plane (0j = (^o) co-polarized (p^'P, s^^s) and 
cross-polarized (p— >s, s— >p) scattering of light incident on the two-dimensional randomly 
rough surface of a dielectric film deposited on the planar surface of silver, whose dielectric 
constant is £3 = —18.28 + 0.4811. The wavelength of the incident light is A = 633 nm, 
the angles of incidence are (60,(^0) = (0.74°, 45°). The dielectric constant of the film is 
£2 = 2. 6896 + 0. Oli, and its mean thickness isd = 478.5 nm. The roughness of the surface 
is characterized by the power spectrum in Eq. l |13| >, with k- = 0.82((b/c), k+ = 1.97(a)/c), 
and its rms height is 5 = A/40 = 15.82 nm. (b) The same as (a) for out-of-plane {(ps = 
(j>0 + 90°) scattering. 



peak lies in the nonradiative region of the q\q2 plane. In addition, among the real satellite peaks 
that should appear in the radiative region, not all may be sufficiently intense to be observable. 



The scattering angles defined by Eq. (14 1 are expected to give the angles at which satellite 
peaks occur in the in-plane co-polarized scattering from the two-dimensional randomly rough 
surface of the film system studied here. 

In the absence of absorption and roughness the {^^(a))} are the solutions of the dispersion 
relation 

a2 (?ll , 0)) = — (e2 [ei jSs (^11 , 0)) + esjSi (^11 , 0))] cot [a2 {q\\,(o)d\ 

± {e| [eift(?|| , (o) + £3)31 {q\\ , 0))] ^ cot^ [a2{q\\ , (o)d] (15a) 

+ 4£i £|£3 j3i {q\\,C0)l5i{q\\,C0)Y^^^ 

in p polarization, and 

a2{q\\,(o) = ^(^ [)3i(^||,«)+)33(^||,a))] cot [a2{q\\,(o)d] 

±{[/3i(^l|,«)+j33(?l|,«)]2cot2 [a2{qii,co)d] (15b) 
+ 4j3i(^l|,«)+i33(?||,«)}'/'^ 



r 9i^/^ 

in s polarization. In these equations j3,(^|| ,Co)= <7| — £i{(o/c) for / = 1,3, while a2(<?|| , O}) 
is defined in Eq. Q. The film structure studied in this paper is found to support two guided 
waves in p-polarization, whose wave numbers are 

^i(a)) = 1.4391(«/c), ^2(«) = l-0119(a)/c), (16a) 

and two guided waves in s polarization, with wave numbers 

= 1.5467(«/c), ^2(«) = l-2432(«/c). (16b) 

These results predict satellite peaks at scattering angles ~ —25.22° and 23.74° in p polar- 
ization and at Qs = —18.13° and 16.65° in s polarization when we are considering in-plane 
scattering, assuming the same angles of incidence as in Fig.[T| These scattering angles are in- 
dicated by vertical dotted lines in Fig. 1(a). The peaks at Qs = —18.13° and 16.65° are seen 
in the s^-s co-polarized scattering contribution to the mean differential reflection coefficient. 

(1 2) 

There is no evidence of satellite peaks at 9s ' — —25.22° and 23.74° in the p^-p co-polarized 
scattering contribution to the mean differential reflection coefficient, presumably because they 
are too weak to be seen. These results disagree with those of |10|, in which no satellite peaks 
were found in the in-plane s— >s scattering contribution to the mean differential reflection coef- 
ficient (although they are present in this contribution when the dielectric film is deposited on a 
planar perfectly conducting surface). However, in 1 10 1 the surface roughness was characterized 
by a Gaussian power spectrum, not the West-O'Donnell power spectrum assumed here. The 
results of earlier calculations p5) of the scattering of p- and s-polarized light from a film with 
a one-dimensional randomly rough surface characterized by a Gaussian power spectrum that 
is deposited on a planar perfectly conducting surface, display satellite peaks more weakly than 
when the roughness is characterized by a West-O'Donnell power spectrum [26]. 

Turning now to the results for out-of-plane scattering presented in Fig. 1(b), we see that an 
enhanced backscattering peak is present in each scattering configuration. It is cut off in each 
configuration. This is an artifact of the present calculation that results from the line defined by 



0i = 00 ± 90° being exactly one grid point away from the backscattering direction. It is im- 
portant to note that in out-of-plane scattering the predominant contribution to the differential 
reflection coefficient is in the cross-polarized part. We see that the satellite peaks are now ob- 
served in the p— >s scattering configuration, meaning that incident p-polarized light excites both 
of the s-polarized guided modes with wave vectors in the = 00 ± 90° directions, which sub- 
sequently interfere to cause satellite peaks in out-of-plane scattering. Hence, the well-known 
"satellite peaks" found in scattering from ID surfaces turn into a kind of "satellite rings" for 
scattering from 2D surfaces, where part of the ring is co-polarized (s— >s in-plane) and part of 
the ring is cross-polarized (p— >s out-of-plane). 

In Fig. 2 we present contour plots of the complete angular distribution of the mean differential 
reflection coefficient for the light scattered incoherently from the film system studied here. The 
material and experimental parameters used in producing these results have the same values used 
in obtaining Fig. 1. 

Light of p polarization (left column) or s polarization (right column) is incident on the struc- 
ture. In Figs. 2(a) and 2(d) all of the scattered light is recorded; in Figs. 2(b) and 2(e) only the 
p-polarized scattered light is recorded; while in Figs. 2(c) and 2(f) only the s-polarized scat- 
tered light is recorded. In Fig. 2(f) we see two regions of high intensity in the in-plane polar- 
ized (s— >s) intensity distribution, centered at radii of approximately 0.29(a)/c) at (j)s = 45°, and 
0.31(a)/c) at 0i = 225°. These are the satellite peaks seen in the plot of {^R<^<^/^Q.s)^^^^^^^ pre- 
sented in Fig. 1(a). No such regions of high intensity are seen in Fig. 2(b) at radii of 0.34(a)/c) 



at 0i — 45° and 0.52(a)/c) at ~ 225°, where satellite peaks are predicted by Eq. ( 14 1 for 
in-plane co-polarized scattering of p-polarized incident light. This result is consistent with the 
absence of satellite peaks in the result for (^dRpp/dQ.s)^^^^^^^ presented in Fig. 1(a). The intensity 
maxima in the out-of-plane cross-polarized (p->s) scattering intensity distribution depicted in 
Fig. 2(c) correspond to the peaks at 9^ « 19° seen in the plot of ^(9/?sp/(9i2v)j^^^j^ presented in 
Fig. 1(b). 

The result that satellite peaks are observed in scattering processes in which the scattered light 
is s polarized, independent of the polarization of the incident light, is an interesting result of the 
present calculations. It may be connected with the fact that s-polarized light is reflected more 
strongly from a dielectric surface than is p-polarized light. 

In Fig. 3 we present results analogous to those presented in Fig. 1, but for angles of incidence 
(0Oj0o) = (5. 19°, 45°). In Fig. 3(a) we present results for the in-plane (0, = <Pq) co-(p— >p, 
s— >-s) and cross-(p— ^s, s— ?>p) polarized scattering, while in Fig. 3(b) we present results for out- 
of-plane (0s = 0oi9O°) co-(p— >p, s— >s) and cross-(p— >^s, s->p) polarized scattering. In the 
results presented in Fig. 3(a) [3(b)] single-scattering processes give no contribution to the mean 
differential reflection coefficient for -46.85° < 0, < 65.57° (-54.59° < 0, < 54.59°). The 
limits of these angular regions are clearly seen in these figures. 

A well-defined enhanced backscattering peak is seen in the results plotted in Figs. 3(a). 



(12) 

Satellite peaks are predicted by Eq. (14 1 to occur (in-plane) at 9,) ' = -31.18° and 19.68° for 



p-polarized incident light, and at 0,'' = —23.20° and 12.30° for s-polarized incident light, 
when the angles of incidence were the same as in Fig. 3. These scattering angles are indi- 
cated by vertical dotted lines in this figure. Peaks at 9s ~ —23.20° and 9s — 12.30° are present 
in the s^-s co-polarized scattering contribution to the mean differential reflection coefficient. 
There is no suggestion of peaks at 9s = —31.8° and 19.68° in the p— >p co-polarized scattering 
contribution to the mean differential reflection coefficient, nor any suggestions of peaks in the 
cross-polarized (p— >s, s— >p) contribution to it. 

The results for out-of-plane scattering presented in Fig. 3(b) show no enhanced backscat- 
tering peaks. The reason for this is simply that since the abscissa points along = 0o ± 90°, 
it does not cut through the backscattering peak, localized at (0^, 0j) — (0o, 0o + 180°). We do 
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Fig. 2. The complete angular distribution of the mean differential reflection coefficient 
{dRg,^ / dSls)mcoh for the light scattered incoherently from the film structure. The material 
and experimental parameters assumed here are those used in obtaining the plots presented 
in Fig. 1. Light of either p (left column) or s (right column) polarization is incident on the 
structure. In (a) and (d) all (diffusely) scattered light is recorded. In (b) and (e) only the 
p-polarized scattered light is recorded, while in (c) and (f) only the s-polarized scattered 
light is recorded. The dark dot in each panel indicates the enhanced backscattering peak. 
Note that the gray scale bar is cut at both ends in order to enhance the sateUite rings. Also 
note that the contribution from single scattering is suppressed, i.e. the differential reflection 

coefficient is artificially set to for qii — kii > fc_ . 




Fig. 3. The same as Fig. 1, but for angles of incidence given by (Oq, (^q) = (5. 19°, 45° 



see some remainders of the satellite ring structure; the low peaks around 9s ~ ±20° are part 
of the rings to the upper left in Fig.|4] As the rings decay in strength away from the direction 
= 00 ± 90°, they are weaker than what is seen in-plane. 

As a necessary, but not sufficient, condition of the validity of our simulation results is energy 
conservation. If all materials in the scattering system are lossless, i.e. Im(e,) = (; = 1,2,3), 
the power of the scattered light has to be equal to the power of the incident light. Under these 
conditions, energy was conserved within 0.03% in our simulations. 
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Fig. 4. The same as Fig. 2, but for angles of incidence given by (%,(^) = (5. 19°, 45°). 
Note that the color bar is cut at both ends in order to enhance the satellite rings. Also 
note that the contribution from single scattering is suppressed, i.e. the differential reflection 

coefficient is artificially set to for qn — kn > A:_ . 



In conclusion, in this paper we have presented a nonperturbative approach to the solution of 
the reduced Rayleigh equation for the scattering of polarized light from a dielectric film with 
a two-dimensional randomly rough surface deposited on a planar metallic surface. We have 
applied this result to calculate the contributions to the mean differential reflection coefficient 
from the in-plane co- and cross-polarized components of the light scattered incoherently, as 
well as from the out-of -plane co-and cross-polarized components of the light scattered incoher- 
ently. The out-of -plane scattering contributions have not been calculated in earlier perturbative 
studies of this problem |10||14) . In addition, we have calculated the full angular distribution of 
the intensity of the scattered light, which has helped to refine the conclusions drawn from the 
calculations of the mean differential reflection coefficient. The main physical result obtained 
in this work is the demonstration that satellite peaks (or rings) can arise in scattering from the 
film structure studied here. This result is in agreement with the results of Soubret et al. fTO) 
but not with those of Kawanishi et al. p4|. A detailed study of the conditions under which 
satellite peaks occur is lacking, but perhaps the approach developed here will be used to deter- 
mine them. The work reported here also opens the door to the possibility of calculating other 
properties of the light scattered from the film system studied here, such as all the elements of 
the associated Mueller matrix, and offers the possibility of designing such structures to possess 
specified scattering properties. These are problems that have to be left to the future. 
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